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1. Introduction 

The single-impurity Anderson model (SIAM) was invented more than forty years ago 
to describe dilute magnetic impurities in metals Ij. In recent years it has become a 
generic model for the physics of strong local electron interactions [2j and is currently 
attracting much interest as a model for studying quantum dots coupled to leads [Sj. 
Its importance has also grown with the advent of the dynamical mean-field theory 
(DMFT) Within the DMFT approach a lattice model such as the Hubbard 

model in the limit of high dimensions is mapped onto an effective SIAM which is 
determined self-consistently. Although the properties of the SIAM are generally well 
understood from numerous studies based on various methods, the situation is not 
fully satisfactory as far as the dynamic properties are concerned, especially the zero- 
temperature high-frequency dynamics. Many-body theories, such as perturbation theory 
or the local moment approach (LMA) provide an explanative (and often accurate) 

picture of quantum impurity dynamics but their accuracy is not known a priori. 
Quantum Monte Carlo (QMC) simulations combined with the Maximum Entropy 
Method jZj give accurate results at finite temperature but extrapolations of the dynamic 
properties to low-temperature are more difficult. The Numerical Renormalisation Group 
(NRG) [HI El EH allows one to determine the low-energy dynamics of quantum impurity 
models almost exactly but this method is less precise at high-energy because of the 
necessary logarithmic discretisation of the host band. 

The NRG has been applied to various quantum impurity problems with great 
success but it has proven more difficult to apply NRG transformations to quantum lattice 
models. These difficulties provided the motivation for the development of the density- 
matrix renormalisation group (DMRG) a decade ago jTTj. DMRG is one of the most 
powerful numerical techniques for studying quantum lattice many-body systems |12j . 
Recently a dynamical DMRG (DDMRG) method has been developed to calculate 
dynamic correlation functions at zero temperature in quantum lattice models [TTH I14j. 
Here we extend this approach to the calculation of the spectral properties in quantum 
impurity models. The method is demonstrated on the flat-band symmetric SIAM. We 
show that DDMRG can provide the spectral density of the SIAM for all frequencies and 
coupling strengths. In particular, the high-energy spectrum can be determined with a 
good resolution. Thus the DDMRG approach is a very useful complement to existing 
techniques for calculating the dynamics of quantum impurity models. 

The Hamiltonian of the SIAM is 

H = E *kfLL + U (n dT - I) (n dl - 1) + £ V k {flX + dtf ka ) , (I ) 

k(J k(T 

where d\ (d a ) creates (annihilates) an electron with spin a =f, J, in a local level (the 
impurity site), hd a = d a d a and fl (f a ) creates (annihilates) an electron with spin a 
in an eigenstate of the (non-interacting) host band with dispersion e k . The sum over 
the index k runs over all states of the host band. The hybridisation between the local 
impurity state and the delocalised band state k is given by the positive couplings 



Density-matrix renormalisation group approach to quantum impurity problems 



3 



Electrons in the local level are subject to a Coulomb repulsion U. (The impurity site 
potential is set by fi = —U/2 as we will discuss the symmetric SIAM only.) 
The impurity one-particle Green's function can be written (77 — ^ + ) 

gm = h 1 A + h — 1 , A , (2) 

\ H - Eq+LO-IT] I \ Eq- H + U + IT] I 

where E is the ground state energy and (...) represents a ground state expectation 
value. In most studies of quantum impurity problems an objective is the computation 
of the impurity spectral density 

D a {u) = --sgn(cu)ImG CT (cu) = A a {uo) + B a {u) (3) 

7T 

with 

AJlu < 0) = lim (dl s ^ A (4) 

B a (u > 0) = lim (4 s ^ Jt \ (5) 

and Ao- (a; > 0) = B a (oj < 0) = 0. The spectral density fulfills the sum rule 

D(uj)Al) = 1. (6) 



The impurity spectral density of the SIAM model depends on the parameter U and 
the hybridisation function 

A(cu) = 7rJ2\V k \ 2 5(u-e k )>0. (7) 

k 

For a symmetric hybridisation function A (a;) = A(— u) the SIAM has a particle-hole 
symmetry. As consequences, the Green's function (J2J) is an odd function of u but the 
spectral density (j3J) is an even function, D a (u) = D a (—oj) [ or A a {uj) = B a {—uS)\, and 
the Fermi energy remains pinned at u = for all U. Furthermore, as the total spin is 
conserved in the SIAM, G a {uj) = G- ff (u). 

In the rest of this paper we will explain how the spectral density of a quantum 
impurity problem such as the symmetric SIAM can be calculated with DDMRG. To 
illustrate the various approaches that we have tested we show results for a well-known 
particular case of the SIAM: the flat-band model Q. In that case, the host band width 
is taken to be much larger than any other energy scales and the hybridisation function is 
assumed to be constant on these scales, A(u) = A . Our goal is then to compute D a {uj) 
in the relevant energy window — W/2 < uj < W/2 with W/2 >U/2,A . It should be 
noted that a Friedel sum rule 

D a (u = 0) = ^- (8) 

7rA 

holds at the Fermi energy u = for all U > in the flat-band symmetric SIAM. In all 
numerical results presented here the energy scale is set by A = 1. 
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2. DMRG method for quantum impurity problems 

Since the invention of DMRG there have been numerous applications of this method 
to systems made of one or more impurities coupled to (generally interacting) one- 
dimensional hosts ^B]- Recently the growing interest for quantum impurity problems 
has spurred the development of DMRG techniques for investigating the more relevant 
problems of impurities coupled to arbitrary non-interacting host bands both in the 
context of quantum dots [THj and of DMFT calculations [13 HE] . 

A direct application of the usual DMRG algorithms to the SIAM Hamiltonian (0) 
is possible but extremely inefficient because it includes "long-range" hopping terms 
(more precisely, an electron can go from one site to any other site in just two hops). A 
better approach is the transformation of the Hamiltonian into a linear chain with 
nearest-neighbor interactions only (as in a NRG calculation [Hj) 

H = V ]T (c\J a + dtc la ) + U (n d] - ~) (n dl - ~) 

+ J]] ajd^Cjfj + Xj (cj a Cj + c^i^Cj^j . (9) 

The new fermion operators Cj a corresponds to electronic states in the host band and are 
related to the original representation by a canonical transformation 

k 

The orthogonal matrix Mjk, the diagonal terms aj and the nearest-neighbor hopping 
terms are calculated with the Lanczos iterative algorithm for tridiagonalising a 
symmetric matrix starting from the initial vector {Mij, = Vk/V} with V 2 = Y^k^k- 
This calculation must be carried out with very high numerical accuracy (quadruple or 
higher precision for floating-point operations) but does not present any other difficulty 
If the original Hamiltonian (0) is particle-hole symmetric, the diagonal terms aj vanish. 

The Hamiltonian © describes an impurity coupled to one end of a one-dimensional 
chain representing the host band states (see figure ^). Here we will use only this 
configuration. However, one can easily imagine different configurations such as an 
impurity site located between two chains as also shown in figure^ In that case, the left 
and right chains could correspond to band states below (e^ < 0) and above (e& > 0) the 
Fermi energy, respectively. Another possibility is that the left and right chains could 
correspond to up-spin and down-spin band states, respectively. The transformation (fTT)|) 
described above can easily be adapted for either cases. It can also be generalised to more 
complicated impurity problems with more than one local level or more than one host 
band. 

The finite-system DMRG algorithm fH\ IT2"] can be used to calculate ground state 
properties of the Hamiltonian (jHJ). In particular, the ground state wavefunction |\l/o) and 
the ground state energy E Q can readily be obtained. To compute dynamic properties 
such as the impurity Green's function we use the dynamical DMRG (Hj. This 
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Figure 1. One-dimensional lattice configurations for applying a DMRG algorithm 
to quantum impurity problems: (a) Impurity site (solid circle) at one end of the chain 
representing the host band (open circles), (b) Impurity site in the middle of two chains. 



approach is based on a variational principle. One can easily show that for rj > and a 
fixed frequency u the minimum of the functional 



(£ + w - H) 2 + rj 1 *\ + 77 (^0 |4| *) + ^ (* |4| ^0 



with respect to all quantum states is 



^(^min) = *C 



—T] 



-dl 



(12) 



[E + u - H\ +r/ z 

The functional minimum is related to the convolution of the spectral density (JSJ) with a 
Lorentz distribution of width r\ by 

W^min) = -nvB*(u). (13) 

A similar result is obtained for the spectral density if one substitutes d a for d' a , —uo 
for uj and A^(uj) for B^(lj) in the above equations. 

The DDMRG method consists essentially in minimising the functional (jllj) 
numerically using the finite-system DMRG algorithm. Thus DDMRG provides the 
spectral densities A%(u) and B^(uj) for a finite broadening 77. The full spectral density 
convolved with the Lorentz distribution 

V 



DKu) 



(14) 



tt[(uj - UJ 1 ) 2 + r/^ 

is given by the sum of A v a {uj) and B%(u). The real part of the Green's function can 
be calculated with no additional computational cost but is generally less accurate. The 
necessary broadening of spectral functions in DDMRG calculations is actually very 
useful for studying continuous spectra or doing a finite-size scaling analysis [Hj . 

The objective of our DMRG simulations is to obtain physical quantities with 
a sufficient accuracy at the lowest possible computational cost. For this purpose, 
the lattice configuration used in the present work makes necessary an adaptation of 
the usual DMRG principles for measurements. In standard DMRG calculations the 
one-dimensional lattice is split in two blocks of sites separated by two intermediate 
sites |11| IT2*] . Calculations are carried out in an effective Hilbert space of dimension 



D ps n fmtTdTii where n is the number of states per site 



; n 



4 counting the spin 



degeneracy in our case) and m,L,R are the numbers of states used to describe the left and 
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right blocks, respectively. Note that rs min(m, A Nl > r ) , where m is the maximum 

number of density-matrix eigenstates kept and Nl and Nr are the numbers of sites 
in the left and right blocks, respectively. The DMRG method errors diminish (often 
exponentially fast) as m is increased while the computational cost increases as m 3 (at 
least in theory). During a DMRG simulation the position where the chain is split is 
moved repeatedly through all sites and the measurement precision for a fixed number 
m varies with that position. Measurements of local quantities such as the density cl a Cj a 
are most accurate if the site j is one of the two intermediate sites. Measurements of 
global quantities, such as the ground state energy or long-range correlations, are most 
accurate if the chain is split in two blocks of equal size because the effective Hilbert 
space dimension D is maximal then. In the present application we have found that 
it is necessary to compromise in order to optimize the ratio between accuracy and 
computational cost in calculations of the spectral density for an impurity located at 
a chain end. Measurements have to be done as close as possible to the impurity but 
far enough from the chain end for the effective Hilbert space to approach its maximal 
dimension D = n 2 m 2 . This implies that measurements are to be done when the left 
block has reached the size of Nl ~ ln(m)/ln(4) sites (assuming the impurity to be at 
the left chain end). 

In practice, we keep enough density-matrix eigenstates to make (D)DMRG 
truncation errors negligible (we have used up to m = 800 states in the present work). 
Thus our DDMRG results for finite systems are numerically exact in the same sense 
as "exact diagonalisation" results are. The main source of errors are finite-size effects, 
which corresponds to the discretisation of the continuous host band in the SIAM and 
are discussed in the next section. 

In figure El we show the spectral density of the flat-band symmetric SIAM calculated 
with DDMRG for {7 = 0. The exact spectral density is a Lorentzian of width A . With 
the broadening (jl4j) it becomes 



On the scale of figure El there is no visible difference between our numerical results 
and this exact result, which demonstrates the accuracy of our method. Note that the 
{7 = limit is a relevant accuracy test. The hybridisation A (a;) with the host band is 
renormalised in the density-matrix renormalisation and thus our method could become 
exceptionally accurate for some specific hybridisations [for instance, A(u) = 0] but 
the local Coulomb interaction U is always treated exactly and thus does not affect the 
method accuracy directly. 

With DDMRG the computational cost (memory and processor time) of computing 
D%(u) for a single frequency to is proportional to the number /V of sites in the system. 
In practice, the total processor time for calculating the full spectrum grows as ~ N 3 
because one usually reduces r\ and increases the frequency resolution as N is increased. 




(15) 
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Figure 2. Spectral density at U = calculated with a constant host band 
discretisation (solid circles) for W = 20A ,-/V = 59, Ae « 0.34A , and ry = 0.5, then 
deconvolved (open circles). Solid lines show the exact results (|15|l without broadening 
and with r\ — 0.5 for comparison. 



3. Discretisation of the host band 

We are interested in the properties of the SIAM with a continuous host band and a 
continuous hybridisation function A (a;) but DMRG calculations can be performed on 
finite lattices only. Therefore, we must discretise the host band and carry out DMRG 
calculations for a finite number N of host band eigenstates corresponding to energies 
€k(k = 1, . . . , N), then extrapolate the results to a continuous host band (N — > oo). 
Choosing an appropriate discretisation of the host band (i.e., selecting the N band 
state energies e^) turns out to be the greatest difficulty in applying DDMRG to the 
SIAM. 

The spectral density (@J) [or similarly (jSJ)] can be written in the Lehmann 
representation 

n 

where \^? n ), n = 0, 1, . . ., design the eigenstates of H contributing to A a {uj) (i.e., with 
non-vanishing matrix elements) and E n their respective energies. On a finite lattice 
this spectrum is always discrete. To determine the exact SIAM spectral density for 
N — > oo it is necessary to use a broadening r\ larger than the space E n+1 — E n between 
two consecutive eigenstates in dense parts of the spectrum (for instance, in continuous 
parts of the spectrum) ^3]. In the flat-band SIAM (and also in self-consistent SIAM 
derived in DMFT calculations for the Hubbard model ^7j) the distribution of excited 
states contributing to the spectral density for finite N is essentially determined by the 



d a 



5(u + E n -E 



(16) 
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distribution of selected energies e^. Therefore, we must use 



T] > E n+1 - E,, 



(17) 



In other words, the host band discretisation directly limits the resolution of spectral 
density calculations. Note that this is a general difficulty for all methods based on a 
discretisation of the SIAM, not only for DDMRG calculations. For instance, a similar 
problem limits the resolution of NRG calculations at high energy. 

In practice, we choose N values e^ in the relevant energy window \u\ < W/2 and 
calculate the corresponding hybridisation terms 



where 5k = (e^-i + £fe)/2. If A (a;) varies slowly for uo pa on the scale Ae& = 
(ejfc+i — ejfc-0/2 then V k 2 pa A(co> = e k )Ae k /ir. For the flat-band SIAM that we are 
considering here this last equation is obviously exact. To preserve the particle- hole 
symmetry in the discrete SIAM we only include pairs of band states with opposite 
energies in the N selected values. For even N this means that there are N/2 different 
pairs (e^ — e^), while for odd N there are (N — l)/2 such pairs and a state at e^ = 0. 
Including the impurity site, the total number of sites in the lattice is thus N + 1. The 
ground state contains an equal number of electrons and has a minimal total spin S = 
for odd N and S = 1/2 for even N. Note that in that last case (S ^ 0) we need to 
compute the spectral density for both a =f and a =[ and take the average, as the 
spin-flip symmetry a <-> —o is broken. It is thus more efficient to work with an odd 
number N of host band states. 

The simplest discretisation scheme consists in choosing N equidistant energies e& 
(Aefc = Ae pa W/N) in the relevant energy window (W/2 > \ej.\). In that case we use 
a constant broadening rj ~ Ae. This approach has been used for the U — results 
shown in figure HJ A constant discretisation is sufficient provided that the spectrum 
does not contain any structure with a width smaller than Ae ex 1/N, which is readily 
achieved in the weak-coupling SIAM. It is well-known [21 EI] that for intermediate to 
strong couplings U the broad spectral feature around u> — shrinks to a sharp peak (the 
so-called Abrikosov-Suhl resonance). The spectral weight is progressively transfered to 
two Hubbard satellites around uj ~ ±U/2 for increasing U. For instance, in figure H3 we 
show the spectral density for the intermediate coupling U = 2.57rA calculated with a 
constant discretisation Ae ~ 0.34A and a broadening 77 = 0.5A . This spectrum agrees 
qualitatively with LMA and NRG results (HI CHI- Quantitatively, however, the DDMRG 
spectrum is obviously inaccurate, especially around to ~ 0. For instance, at the Fermi 
level the Friedel sum rule (JHJ is clearly not fulfilled. As for the U = case this apparent 
discrepancy is due to the broadening rj. 

To obtain better results one can perform calculations with a higher resolution 
(which means smaller 77 pa Ae ~ 1/N and thus larger lattice size N) and possibly 
extrapolate to the limits 77 — > and N — > 00. In figure we show D%(ui = 0) calculated 
with DDMRG at U — 2.57rAo for several values of rj. For comparison we also show the 




(18) 
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Figure 3. Spectral density for U = 2.57rA and W = 24A . Calculated for a 
constant discretisation Ae rj 0.34A o and i] — 0.5A (dashed line), calculated using 
variable discretisation 0.2 > Ae/Ao > 0.0059 and broadening 0.25 > t]/Aq > 0.01 
(solid line), and calculated using a variable discretisation 0.76 > Ae/Ao > 0.16 and 
a constant broadening r\ = 0.25Ao then deconvolved (circles). Inset: expanded view 
around the Fermi level u> = 0. 

exact scaling for the non- interacting (U = 0) case. Clearly, D%(u = 0) tends to the 
exact result as rj — > but the convergence is slow and will become worse for sharper 
resonances (larger U). For instance, for U = 87rA we can already clearly observe a 
very sharp Abrikosov-Suhl resonance using a resolution of Ae m r] — 10 _4 Ao but the 
height of the peak is only D^(uj = 0) ~ 0.06/7rA . Actually, one can guess that a 
broadening 77 smaller than the Kondo scale cx exp(— ttU/8Aq) [21 E] is required to obtain 
the low-energy spectral density with a good resolution and accuracy. Therefore, the 
required system size N ~ rj' 1 oc exp(7r?7/8Ao) increases exponentially with U if one 
uses a constant discretisation. This approach is clearly not practicable. 

As an interest of many quantum impurity problems is the investigation of sharp 
resonances in the spectral density and a analogous sharp quasi-particle peak is an 
essential feature of DMFT calculations, it is desirable to achieve a high accuracy and 
resolution for these spectral features too. A better approach is the use of a variable 
discretisation. For instance, one can select two sets of equidistant energies in the 
host band with a energy resolution Ae for high energies W/2 > je^l > W*/2 and with 
another energy resolution Ae' < Ae for |efc| < W*/2 (see figure EJ. In that case the 
broadening rj also depends on u. Typically, we use rj w Ae for W/2 > \u\ > W* /2 
and 7] ~ Ae' for < W*/2. This allows us to obtain the spectral density around 
the Abrikosov-Suhl resonance with a much higher resolution for a given number of host 
band sites N. It is also possible to use more than two different energy resolutions, to 
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Figure 4. Fermi level spectral density ttAoD2-(lo = 0) as a function of the broadening 
r\. Exact result for U = (dashed line) and DDMRG results for U = (squares) and 
[/ = 2.57rA (circles). 



use a higher energy resolution for the high energy sector than for the low energy sector, 
or even to select a higher resolution at intermediate energies (see figure EJ). Moreover, it 
is possible to combine the results of the high-resolution sectors obtained with different 
variable discretisations in a single high resolution spectrum. 

This approach allows us to achieve a resolution which is significantly better than 
with the constant discretisation for the same number of host band states N. For 
instance, figure El shows the spectral density obtained with variable discretisation Ae 
and broadening rj ranging from Ae = 0.2A (rj = 0.25A ) for the Hubbard satellites to 
Ae = 0.0059A ( rj = 0.01A ) for u — > 0. This approach clearly gives much better results 
for the sharp Abrikosov-Suhl resonance than the constant discretisation approach. The 
Friedel sum rule (JSJ) at the Fermi level is fulfilled within 3 %. It should be noted that 
this approach is not without problem. The most obvious drawback is that a variable 
broadening rj = i](uj) breaks the sum rule © and can significantly change the spectrum 
shape if it varies faster than the bare spectrum as a function of uj. (This could be the 
case for the Abrikosov-Suhl resonance shown as a solid line in the inset of figure 3.) 

To take advantage of the higher resolution provided by a variable discretisation 
while keeping the benefit of a constant broadening rj we have devise the following scheme. 
The relevant energy window is split into several intervals. For each interval, a DDMRG 
calculation is done using a variable discretisation with a high resolution in that interval 
and a lower resolution outside. The spectral function is calculated for frequencies in 
that interval. Then the results obtained in each calculation are combined into a full 
spectrum. The computational effort is equal to that of a single calculation with a 
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Figure 5. Various host band discretisation schemes: (a) constant, (b) variable with 
a higher resolution at low energy, (c) variable with a higher resolution at intermediate 
energy, and (d) logarithmic. 



constant discretisation for the same number of host band states, but the resolution is 
higher. In figure El we show the spectral density calculated using this scheme with a 
high resolution of A t = 0.16 — 0.20A inside each interval and a lower resolution of 
A e = 0.71 — 0.76A outside. The spectrum has been calculated with r) = 0.25A then 
deconvolved (see next section). This spectrum is clearly sharper than the one obtained 
with a constant discretisation and, by comparison with related NRG and LMA results, 
it seems also more accurate. The improvement is due both to the lower broadening r\ 
and the deconvolution. Nevertheless, the variable discretisation scheme with varying 
broadening rj is more accurate for low frequency and gives a much sharper Abrikosov- 
Suhl resonance. 

As in NRG calculations [BJ a logarithmic discretisation of the host band tk = 
{W/2)A- k (with A > 1 and k = 1,2,..., N/2) is possible. Ground state DMRG 
calculations can easily been performed although the computational effort is significantly 
greater than with NRG. Calculating the spectral density with DDMRG is more 
problematic. The main problem is that we have not found any satisfactory method 
to broaden the spectrum. In NRG calculations [32] the spectrum is convolved with a 
function which vanishes for uj — > and broadens the spectrum on a logarithmic scale 
rj oc uj ~ Aefc. In DDMRG calculations we always get the spectrum convolved with a 
Lorentzian (|14j). To imitate the logarithmic broadening we have tried a "Lorentzian" 
broadening (JT3j) with a variable rj = T)(oj') = Xuo' . Note that with the DDMRG method 
we can change rj as a function of uj but not of uj' in (JT3J). Fortunately, if G a (uj) is the 
impurity Green's function obtained with DDMRG for T](u) = Xuj, the desired Green's 
function G a (u) with t](uj') = Xuj' can be obtained through 



Note that the condition ()17|) implies A > 1 — A -1 . We have tested this approach on 
the flat-band SIAM for various values of U. As expected, the high-energy spectrum 
is widely broadened and qualitatively similar to corresponding NRG results [TU] . In 
the low-energy spectrum, however, this scheme clearly does not work because D a {u) 
seems to diverge for uj — > or at least the Fermi level spectral density -0^(0) is much 
larger than the exact value l/(7rA ). We think that this failure is due to the broadening 




(19) 
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with the Lorentzian function 1)140 . which does not vanish for uo — > contrary to the 
function used in NRG calculations [10J. Therefore, a logarithmic discretisation does not 
seem to be useful for DDMRG calculations. It should be noted that with the Lanczos- 
DMRG 20 or the correction-vector DMRG [2T] it should be possible to calculate the 
spectral density in the Lehmann representation ()TH)l (i.e., without broadening) and to 
broaden it on a logarithmic scale as in NRG calculations. This approach could give 
better results for the low-energy spectrum than the DDMRG approach discussed here 
but has yet to be tested. 

In summary, the host band discretisation determines the resolution of DDMRG 
calculations for the impurity spectral density. No single discretisation scheme is 
appropriate for all cases. Nevertheless, a combination of different schemes can be 
used for the various features of the same spectrum. Therefore, DDMRG allows one 
to calculate a spectral density with high resolution for all frequencies provided the host 
band discretisation and the broadening are adapted to the specific problem and its 
spectral features. 

4. Deconvolution of DDMRG spectra 

An approach for obtaining sharper spectra is the deconvolution of the DDMRG data. 
In theory, a deconvolution means solving ()14)) for D(u) using the DDMRG data in the 
left-hand side. However, this is typically a ill-conditioned inversion problem ^H] which 
generally cannot be solved numerically. Moreover, if it was possible to do this calculation 
exactly, one would obtain the discrete spectral density of the SIAM on a finite lattice 
of N + 1 sites. 

Nevertheless, the broadened spectral density of the SIAM on a infinite lattice 
(N — > oo) is usually almost identical to the spectral density of the discretised SIAM 
(N < oo). For instance, for U — there is no visible difference in figure 121 between 
the DDMRG results for N = 59 and the exact result for N — > oo. [A necessary 
condition for a large enough n is (117)) .] Therefore, one can make the approximation that 
DDMRG data for D%(lj) describe the broadened spectral density for iV — > oo and solve 
(fT4")l approximately under the condition that D a (u) is the exact spectral density of the 
SIAM. For instance, one can require that D a (u) is a continuous and relatively smooth 
function. 

In practice, we calculate D a {uo) using the algorithm presented in or using 
standard linear regularisation methods for the inverse problem ^H]- Other possibilities 
include using Fourier transformations (THj or the Maximum Entropy Method as in QMC 
simulations 0. In figure |2] one sees that the deconvolved DDMRG spectral density 
agrees perfectly with the exact result for U = 0. Similarly, a deconvolution gives 
excellent results for finite but weak coupling U. For U = 2.57rA the deconvolved 
spectral density in figure 01 shows sharper Abrikosov-Suhl resonance and Hubbard 
satellites than the original DDMRG data for rj = 0.25A (not shown). However, the 
Friedel sum rule (JBJ) at the Fermi level is far from being fulfilled (the relative error is 
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Figure 6. Upper Hubbard satellite for U = 8ttA . Calculated with DDMRG for 
W = 40 A and a constant discretisation Ae = 1.29A and 77 = A (dashed line) then 
deconvolved (circles). The solid line shows the LMA result [H] for U :» Ao (see text). 

about 10 %), which indicates that the Fermi level resonance is still too broad despite 
the deconvolution. In general, we have found that the deconvolution works very well 
for the Hubbard satellites and other broad structures (broader than 77) but bring only 
a partial improvement for sharp peaks. 

Deconvolution methods become rapidly instable if a variable broadening n is used 
in Q14|) . Therefore, we apply them to spectra calculated with a constant n only. 
Nevertheless, if the spectrum is made of well separated structures (as the Abrikosov- 
Suhl peak and the Hubbard satellites are for strong coupling U ^> A ), it is possible 
to deconvolve each structure separately. In that case we can use a different n for each 
structure. In figure El we show the upper Hubbard satellite calculated with DDMRG 
for U = 8ttAq then deconvolved. Our numerical result agrees remarkably well with 
the LMA prediction for this spectral structure in the strong coupling regime [/ > A : 
a Lorentz distribution of width 2Aq and total weight \ centered at u c = U/2. Note 
that in figure |B1 we have chosen uj c = 13.7Ao to fit our data rather than the value 
u c = U/2 = 47rA ~ 12.6A . This shift of 1.1A is due to corrections to uj c of the 
order of A /U |Hj and the finite discretisation Ae = 1.29A used in our calculation. The 
deconvolved DDMRG spectral density appears to contain more spectral weight than the 
LMA result because we have used a rather narrow energy window < W/2 = 20Ao 
in the DDMRG calculation and thus the spectral weight which normally lies in the 
high energy tail (|u;| > W/2) of the Lorentz distribution has been shifted to low energy 
(M < W/2). 
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5. Conclusion 

We have extended the DDMRG method ^3] to the investigation of dynamic properties 
in quantum impurity problems. The method has been demonstrated on the flat-band 
symmetric SIAM. We have obtained accurate results for the impurity spectral density 
for all frequencies and coupling strengths. The main difficulty of this approach is the 
choice of an appropriate discretisation of the host band to resolve sharp spectral features 
such as the Abrikosov-Suhl resonance in the strong-coupling regime of the SIAM. 

Our method can readily be used to study quantum lattice many-body systems in 
the framework of the DMFT. We have already applied this DDMRG-DMFT technique 
to the Hubbard model in infinite dimensions. For the metallic phase at weak coupling 
U our numerical results have a better resolution than those obtained with NRG and 
are in excellent agreement with weak-coupling perturbation theory In the Mott 

insulating phase DDMRG-DMFT calculations are relatively easy to perform because 
there is no sharp structure to resolve around uo = 0. Thus we can obtain very accurate 
results for the Hubbard bands in the spectral density and determined the one-particle 
gap. Our numerical results agree remarkably well with strong-coupling perturbation 
theory [22]- We are currently investigating the intermediate-coupling regime, where the 
Mott metal-insulator transition occurs [22]) using DDMRG-DMFT. 

The DDMRG method presented here can be generalised and applied to a great 
variety of quantum impurity problems. First, although we have discussed only the 
symmetric flat-band SIAM, our method can readily be applied to an asymmetric 
SIAM [6] or more complicated hybridisation functions. Second, the method can be used 
for other Hamiltonians than the SIAM (for instance, the Kondo model) provided that the 
transformation to a one-dimensional lattice configuration with short-range interactions 
only is possible. Third, contrary to NRG our method can resolve sharp spectral 
structures at any frequency and thus could be very useful for studying the dynamics 
of systems with resonances at finite frequency, such as the Kondo model in magnetic 
field (23] . Fourth, DDMRG can be applied to models with other degrees of freedom than 
fermions. For instance, there are very efficient density-matrix renormalisation methods 
to treat bosons [23]. One can easily combine them with DDMRG to investigate a 
quantum degree of freedom coupled to a dissipative environment such as a phonon 
bath [2B]. Finally, our method can be generalised to problems with more than one 
impurity or more than one host band. If the impurity sites are arranged in a one- 
dimensional lattice and the number of host bands do not exceed two, the computational 
effort is probably comparable to that required for the investigation of the SIAM 
presented here. Therefore, one should be able to investigate the dynamics of a two- 
channel Kondo problem or of quantum dots and wires coupled to (up to) two leads [3] 
without difficulty If the impurity sites form other structures such as a two-dimensional 
cluster or the number of host bands is greater than two (which occurs for multi- 
channel Kondo problems or in the dynamical cluster approximation 27 ]), DDMRG 
calculations of dynamic properties are also possible but the computational effort could be 
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substantially greater. In summary, DDMRG provides a powerful and versatile approach 
for investigating the dynamics of quantum impurity problems. 
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